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ABSTRACT 


« 


A  number  of  experimental  studies  are  being  conducted  to 
validate  theoretical  methods  of  calculating  fallout  protection 
afforded  by  structures.  As  a  part  of  this  program,  exposure 
rates  in  an  open  basement  were  measured  at  the  US  Army  Nuclear 
Defense  Laboratory  (USANDL).  This  report  presents  calculated 
exposure  rates  for  forty-eight  detectors  located  within  the 
open  basement.  These  results  were  obtained  using  UNC-SAM-2, 
a  Monte  Carlo  radiation  transport  digital  computer  code.  The 
calculated  results  are  compared  with  experimental  measurements 
and  also  with  comparable  adjoint  (Monte  Carlo)  results. 


FOREWORD 

It  is  a  pleasure  to  acknowledge  the  contributions  of 
Dr.  L.  M.  Petrie  for  his  initial  work  with  the  UNC-SAM-2  code 
and  Mr.  A.  T.  Futterer  for  his  assistance  and  continued  interest 
is  this  work. 
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CALCUIATIONS  WITH  UNC-SAM-2  OF 
EXPOSURE  RATES  MEASURED  IN  AN  OPEN  BASEMENT 

1.  BACKGROUND 

The  US  Army  Nuclear  Defense  Laboratory  (USANDL)  haB 
conducted  experiments  to  determine  the  shielding  afforded  by 
open  and  concrete-covered  basements  (References  1  and  2).  The 
theoretical  methods  set  forth  by  Spencer  (Reference  3)  and  the 
engineering  methods  presented  by  Eisenhauer  (References  L  and  5/ 
were  used  by  USANDL  to  calculate  protection  factors  for  both 
these  experiments. 

This  report  presents  an  effort  to  calculate  theoret¬ 
ically  the  exposure  rates  for  the  open  basement  experiment 
(Reference  l)  using  Monte  Carlo  techniques. 

This  approach  is  necessary  in  order  to  isolate  certain 
factors  of  fallout  protection  which  are  not  amenable  to  direct 
experimental  verification. 

The  Monte  Carlo  method  is  best  suited  for  the  calcula¬ 
tion  of  fallout  protection  afforded  by  structures  for  two  reasons. 
First,  more  complex  geometries  can  be  represented  via  Monte 
Carlo  than  by  other  methods.  Second,  factors  associated 
with  the  engineering  methods,  such  as  exposure  contributions 
due  to  skyshine,  in  and  down,  etc.,  are  easily  separated  in 
Monte  Carlo  calculations  by  placing  infinite  absorption  cross 
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sections  in  certain  volumes  and  allowing  only  certain  desired 
contributions.  The  major  difficulty  associated  with  the  fallout 
problem  is  the  very  large  magnitude  of  the  ratio  of  source  to 
detection  area.  One  suggested  solution  to  this  problem  is  the 
use  of  an  adjoint  Monte  Carlo  method  in  which  particles  are 
traced  '.w  kward  from  detector  to  source. 

The  present  report  treats  preliminary  work  done  in 
the  Theoretical  Physics  Branch  of  USANDL  and  is  concerned 
exclusively  wiu  the  open-basement  program.  A  Monte  Carlo 
code,  UNC-SAM-2  (Reference  6)  was  used  to  calculate  exposure 
rates  for  forty-eight  detectors  in  the  open  basement.  These 
exposure  rates  are  compared  with  experimental  results  and  with 
adjoint  Monte  Carlo  calculations  from  the  GADJET  code  (Refer¬ 
ence  7). 

Preliminary  calculations  are  beinr  performed  for  the 
closed  basement  problem,  and  the  results  will  soon  be  forth¬ 
coming.  Further  work  will  be  proposed  to  make  the  open  and 
closed  basement  calculations  for  bettei  accuracies  and  also 
to  isolate  the  sources  of  the  various  components  of  the  exposure 
rate  for  validation  or  modification  of  the  engineering  methods. 

2.  EXPERIMENT 

Exposure  rates  were  measured  at  USANDL  (Reference  l) 
in  the  open  basement  shown  in  Figure  2.1.  Beginning  1  foot 
belcw  the  top  of  the  structure,  detectors  were  spaced  vertically 
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Figure  2.1  Layout  of  the  basement  and  fallout  field. 

1  foot  apart  in  banks  of  six  (Figures  2.2  and  2.3).  Forty- 
eight  detectors  were  located  in  banks  at  positions  A,  B,  C,  D, 

E  ,  as  ,  bs  ,  and  es  as  indicated  in  Figure  2.b.  A  s°Co  source 
■was  circulated  within  the  tubing  seen  the  background  of 
Figure  2.1.  The  tubing  was  laid  in  a  semicircular  configura¬ 
tion  on  the  180°  flat  field  adjacent  to  the  basement.  The 
experiment  was  run  for  several  annular  radii;  however,  the 
present  calculations  are  for  a  source  radius  of  60  feet  (Refer¬ 
ences  1  and  3).  Also,  for  calculation  purposes,  it  is  assumed 
that  the  source  is  uniformly  distributed  (l  Ci 
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Figure  2.4  Top  view  of  the  basement  showing  primary  and  image 
detector  locations. 

This  experiment  was  chosen  as  a  test  for  the  calcula¬ 
tions!  method  because  of  its  simple  geometry,  accessible  exper¬ 
imental  results  and  available  adjoint  (GADJET)  calculated 
results  (Reference  7)*  Little  difficulty  would  be  encountered 
in  performijig  similar  calculations  for  other  structures  if  the 
geometries  were  relatively  simple  and  cross  sections  were 
available . 
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3.  CALCULATIONAL  METHOD 


The  Stochastic  Approximation  Method  (SAM)  developed 
by  United  Nuclear  Corporation  is  well  documented  (Refeience  6) 
and  will  not  be  discussed  here  in  detail.  Only  the  SAM 
information  required  to  calculate  the  open  basement  results 
will  be  included  in  this  report. 

The  salient  specifications  for  UNC-SAM-2  are : 

(1)  Source:  Isotropic  60Co  emitters  were  uniformly 

distributed  in  a  thin  air  layer  (10~6  cm) 
at  the  air-ground  interface.  The  source 
energy  peaks  were  1.17  and  1.33  MeV. 

(2)  Transport:  The  interactions  allowed  were  Compton 

scattering  and  absorption  by  photoelectric 
and  pair  production  processes. 

(3)  Detection:  Point  detectors  and  the  Flux  at  a  Point 

(FAP)  technique  were  used.  The  lower 
energy  cutoff  of  these  detectors  was 
20  keV. 

The  geometry  setup  for  the  basement  is  a  straight¬ 
forward  assembly  of  thirty- three  geometric  regions.  For 
example ,  each  wall  of  the  basement  was  represented  by  a  con¬ 
crete  box  with  the  appropriate  dimensions ;  the  atmosphere  con¬ 
sisted  of  several  large  boxes,  etc.  Point  detectors  were 
placed  at  the  appropriate  coordinates  within  the  basement 
region. 


The  fifteen  energy  groups  shown  in  Table  3»1  were 
used  throughout  the  computations. 
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TABLE  3.1  ENERGY  GROUP  STRUCTURE 


Group 

Energy  (MeV ) 

Group 

Energy  (MeV) 

T 

A. 

1.35  - 

1.25 

9 

o 

i 

.50 

2 

1.25  - 

1.15 

10 

.50  - 

.40 

3 

1.15  - 

1.05 

11 

.40  - 

.30 

4 

1.05  - 

1.00 

12 

.30  - 

.20 

5 

1.00  - 

.90 

13 

.20  - 

.10 

6 

.90  - 

.80 

14 

.10  - 

.05 

7 

.80  - 

.70 

15 

.05  - 

.02 

8 

.70  - 

.60 

Both  regional  and  angular  weighting  schemes  were  used 
in  an  effort  tc  reduce  the  variance.  These  were  established  in 
such  a  way  as  to  give  photons  moving  towards  the  open  basement 
more  importance  then  the  photons  moving  away  from  the  basement. 
However,  in  this  initial  calculation,  no  effort  was  made  to 
optimize  the  weighting  scheme. 

As  indicated,  Flux  At  a  Point  (FAP)  technique  (Refer¬ 
ence  9)  was  used  for  the  energy  dependent  detector  flux  calcula¬ 
tion.  Flux  values  were  also  determined  for  certain  physical 
regions . 

Subroutine  SCR  was  written  to  calculate  flux  to 
exposure  as  follows : 
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D(E)  =  &H(E)Eh.(e)y(E)aE  (R  h-1  ) 


“total  -  I  “<E> 
all  A E 


where : 

3  =  constant 

H(E)  =  detector  response  function 

E  -  average  energy  per  energy  group,  (MeV) 

^,(E)  =  air  absorption  cross  section  (cm2) 

(E )  =  particle  track  length  per  cm3  per  E  per 

initial  particle 

AE  =  group  width  (MeV). 

H(E)  for  the  Victoreen  detectors  is  shown  in 
Figure  3*1  (Reference  8). 

To  obtain  the  total  exposure  from  a  360°  fallout 
field,  exposure  rates  were  calculated  for  three  detector  banks 
located  at  positions  a«  ,  t.  ,  and  ea  (Figure  2.4).  The 
exposure  rates  calculated  for  the  180°  fallout  field  at  posi¬ 
tions  C  and  D  are  doubled  to  obtain  the  exposure  rates  for  the 
360°  field.  However,  due  to  nonsymmetry,  the  exposure  rates 
for  a  360°  fallout  field  at  positions  A  ,  B  ,  and  E  ,  are 
obtained  by  adding  doses  calculated  for  the  l80°  field  as 
follows : 
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Figure  3.1  Energy  response  curve  for  Victoreen  H(E)  detectors 
(Reference  8). 

A  (360° )  =  Ai  (180° )  +  as (180° ) 

B  (360° )  =  Bi (180° )  +  b2 (180° ) 

E  (360° )  =  Ei (180° )  +  es (180° ) 

The  GADJET  method  as  well  as  the  GADJET  results  pre¬ 
sented  here  are  well  documented  (Reference  7);  suffice  it  to 
say  that  the  GADJET  code  is  based  on  solving  the  adjoint 
Boltzmann  transport  equation.  The  code  traces  particles  from 
the  detector  to  the  source.  All  photons  fo:  which  calculations 
are  performed  contribute  to  the  final  result. 
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4.  RESULTS  AND  COMPARISON 

The  experimental,  calculated,  and  GADJET  exposure 
rates  are  given  in  Table  4.1  for  each  of  the  thirty  detectors 
located  in  banks  of  six  each  a  positions  A^  ,  Bi ,  C,  D,  and  Ei 
in  Figure  2.4.  These  results  are  plotted  versus  height  above 
the  basement  floor  in  Figures  4.1  through  4.5  for  these  posi¬ 
tions  . 

Except  for  the  detectors  near  the  top  surface,  the 
UNC-SAM-2  results  sre  in  agreement  with  the  experimental 
results.  Note  that  experimentally  determined  exposure  rates 
for  the  detectors  near  the  top  of  the  basement  are  consider¬ 
ably  higher  than  those  calculated  by  either  UNC-SAM  or  GADJET. 

The  calculated  standard  deviation  for  the  SAM  results 
is  approximately  30  percent  for  a  given  detector.  Only  region- 
angle  weighting  schemes  were  used  in  this  initial  investiga¬ 
tion,  and  no  optimum  biasing  schemes  were  used. 

This  work  has  shown  that  UNC-SAM-2  can  be  used 
effectively  for  fallout  shielding  calculations.  The  statistical 
accuracy  could  be  improved  by  ircreasing  the  number  of  histories 
or  with  the  use  of  more  effective  biasing.  Further  work  will 
be  done  to  isolate  the  sources  of  the  various  components  of  the 
exposure  rate  for  validation  or  modification  of  the  engineering 
method . 
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TABLE  4.1  COMPARISON  OF  EXPOSURE  RATES  OF  DETECTORS  AT  VARIOUS  LOCATIONS  IN  AN  OPEN 
CONCRETE  BASEMENT 


5 

S 


Note i  Exposure  rates  are  in  mR  h 

For  detector  locations,  see  Figures  2.3  and  2.^. 


Detector 

Height 

Above  FI'  or 
(ft) 

Exposure 

Rates 

Experi¬ 

mental 

Calcu¬ 

lated 

(SAM) 

Adjoint 
(GAB JET) 

Experi¬ 

mental 

Calcu¬ 

lated 

(SAM) 

Adjoint 
(GAD JET) 

Detector  Location  A 

Detector  Location  B 

1 

666 

651 

464 

880 

759 

£38 

2 

781 

708 

304 

1050 

909 

892 

3 

826 

1244 

537 

1240 

1539 

1013 

4 

1010 

1121 

890 

1450 

1375 

1419 

5 

1320 

1949 

1470 

i960 

2260 

2057 

6 

2340 

2208 

1665 

3360 

2995 

2345 

Detector  Location  C 

Detector  Location  D 

1 

950 

682 

879 

754 

780 

509 

2 

1090 

1560 

949 

920 

1092 

500 

3 

1320 

1717 

1311 

992 

1170 

764 

4 

1560 

2030 

1219 

1170 

IO92 

743 

5 

2000 

2557 

1785 

1700 

1873 

1368 

6 

3900 

3016 

2488 

2820 

2341 

2175 

Detector  Location  E 

1 

864 

760 

777 

2 

976 

846 

767 

3 

1160 

1215 

902 

4 

1310 

1216 

1448 

5 

1750 

1661 

1323 

6 

3700 

1981 

1482 

S 


Figure  U.l  Exposure  rates  for  detectors  at  position  A  in  an 
open  basement. 


Figure  b.2  Exposure  rates  for  detectors  at  position  B  in  an 
open  basement. 
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